Clustering of Ions at Atomic Dimensions in Quantum Plasmas 
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By means of particle simulations of the equations of motion for ions interacting among themselves 
under the influence of newly discovered Shukla-Eliasson attractive force (SEAF) in a dense quantum 
plasma, we demonstrate that the SEAF can bring ions closer at atomic dimensions. We present 
simulation results of the dynamics of an ensemble of ions in the presence of the SEAF without and 
with confining external potentials and collisions between the ions and degenerate electrons. Our 
particle simulations reveal that under the SEAF, ions attract each other, come closer and form ionic 
clusters in the bath of degenerate electrons that shield the ions. Furthermore, an external confining 
potential produces robust ion clusters that can have cigar-like and ball-like shapes, which remain 
stable when the confining potential is removed. The stability of ion clusters is discussed. Our results 
may have applications to solid density plasmas (density exceeding 10 23 per cubic centimeters), where 
the electrons will be degenerate and quantum forces due to the electron recoil effect caused by the 
overlapping of electron wave functions and electron tunneling through the Bohm potential, electron- 
exchange and electron-exchange and electron correlations associated with electron-1/2 spin effect, 
and the quantum statistical pressure of the degenerate electrons play a decisive role. 
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During the early thirties, there were several discover- 
ies related to non-Coulombic shielded potential distri- 
butions that exhibit the role of collective interactions 



between electrons and ions in electro-chemistry (Debye 
& Hiickel 1923[ ) (viz. electrolytes and co l loidal sus 



tures under the influence of Coulombic, DH, TF and 
Yukawa repulsive forces, one must confine the like- 
charged particles in an external potential, so as to bring 
them to a minimum energy state. Examples include the 
Wigner crystals (Wigner 1934) composed of an ensem- 



pensions), solid state (Fermi 1927 Thomas 19271 and ble of electrons on the surface of liquid helium, ion crys- 



gaseous ( |Langmuir| 1929) plasmas, and between neutrons tals in laser cooled Paul (Drewsen et al. 1998) and Pen 



and protons in elementary particle physics (Yukawa 
1935). The screened non-Coulombic potentials, which 
were obtained by using linearized theory based on the as- 
sumption that the potential energy between the particles 
is much smaller than the particle kinetic energy, are now 
known as the Debye-Hiickel (DH), Thomas- Fermi (TF), 
and Yukawa potentials in the context of electro-chemistry 
and plasma physics, condensed matter physics, and nu- 
clear physics, respectively. The DH, TF, and Yukawa 
potentials describe short-range (of the order of the DH 
radius, the TF radius, and the Yukawa radius, which 
are fixed by the size of a shielded cloud) repulsive in- 
teractions between two particles that have the same po- 
larity. The DH theory has also been extended to dusty 
plasma physics where charged dust particles are shielded 
by non-degenerate electrons and ions. The DH, TF, 
and Yukawa interaction potentials, which significantly 
deviate from the long-range Coulomb interaction poten- 
tial, have important applications to the understanding of 
phase transitions ( Avinash 2007 Klumov 2010 Kremer 



et al. 1986 ) in different areas of physical sciences. 



In order for charged particles to form ordered struc- 
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ning (electromagnetic) traps (Tan et al. 1995), charged 



dust particle crystals (Wuerker et al. 19591, which were 



formed when like-charged dust particles were kept to- 
gether via external confining potentials despite short- 
range Coulombic or shielded Coulombic repulsive forces 
between charged dust particles. In fact, both electron 
and ion crystals, as well as crystals of colloidal suspen- 
sions and oil droplets have been observed experimentally 



under different physical circumstances (Berg & Gauk- 


lerj [19691 |Deshpande & Bockrath[ |2008 |Drewsen et al. 


1998,|Grimes & Adams[ 1979, M0lhave & Drewsen[ 


2000 


Robertson & Younger 1999 Staanum et al. 2010[| r 


"an et 


al. 


1995 Winter & Ortjohann 1991). Moreover, an en- 



semble of strongly correlated micron-sized negative dust 
particles form dust Coulomb crystals (Barkan & Mer- 
Imol[T995l|Chu fc I[ [T9941 |Fortov et ai\ |19971|Hayashi fc 
Tachibanaj |1994[ |Mohideen et al.\ |1998[ |Thomas et al.\ 
1994) when they were confined by the sheath parabolic 



potential in low-temperature laboratory dusty plasma 



disch arges (Shukla & Eliasson 2009 Shukla & Mamun 
2002). The condensation of charged dust particles oc- 



curs since the dusty plasma Yd (the ratio between the 
Coulomb energy between highly charged dust grains and 
the average dust particle kinetic energy) becomes rela- 
tively large due to the high dust charge state and low 
dust temperature. The attraction between like-charged 
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dust particles forming dust Coulomb crystals may also be 
attributed to attractive forces arising from overlapping of 
the dusty plasma Debye spheres ( Resendes et al. [2998 ) , 
ion focusing and wakefield effectsQNambu et al.\ 1995 



Shukla & Rao 1996J [Vladimirov fc NambuU1995 | ), and 



dust dipole-dipole interactions (Mohideen et al. 1998 



Shukla & Eliassonl 2009 Shukla & Mamunl 120021). The 



Cooper's pairing of charged dust particles, which are 
glued by ions, led to the discovery of a soft- condensed 
matter of dust particle crystals in low-density and low- 
temperature classical plasma with Maxwell-Boltzmann 
distributions for electrons and ions. It turns out that 
several milestones were reached in the areas of ordered 
crystalline structures composed of charged particles (e.g. 
an ensemble of in low-temperature physical systems elec- 
trons, ions, as well as charged colloidal and dust particles) 
in physical systems, which share some common physics. 

However, solid density plasmas are of fundamental im- 
portance for industrial applications (e.g. semiconductors, 
nano-diodes and metallic nanostructures for thin films), 
for inertial confinement fusion (ICF) schemes that uti- 
lizes high density compressed (HDC) plasmas, as well 
as for planetary systems [e.g. the core of Jupiter (For- 



tov 



e.g 



20091] and superdense astrophysical objects 
the cores of white dwarf stars, warm dense matter). 
In dense plasmas, one has to account for degeneracy 
( Chandrasekhar| 1931 |1939[ ) of electrons which obey 



the Fermi-Dirac distribution function. Correspondingly, 
quantum mechanical effects play a vital role since in such 
dense plasmas the Wigner-Seitz radius d = (3/47rno) 1,/3 
is comparable to the thermal de Broglie wavelength 
As = h/mVr (which is a measure of the extent of elec- 
tron wave functions), where h is Planck's constant di- 
vided by 27T, m the electron mass, V? = y/kgT/m 
the electron thermal speed due to random electron mo- 
tions, and ks the Boltzmann constant. Also, in dense 
plasmas with degenerate electrons, Xb turns out to be 
much smaller than the Landau length = e 2 /fcsT, 
which can be conveniently expressed as ksT -C e 2 /as, 
where as = Ti 2 /me 2 is the Bohr radius of a hydrogen 
atom. The electron degeneracy effects at nanoscales in 
dense plasmas can thus be captured through the consid- 
eration of the Fermi-Dirac statistics for electrons with 
spin-1/2 (Fermions), and overlapping of electron wave- 
functions due to Heisenberg's uncertainty principle and 
Pauli's exclusion principle, as well as electron-exchange 
and electron-correlations. Hence, there are quantum 



forces (Brodin et al.. 2008 Chandrasekhar 1931 


1939 


Crouseilles et al.\ 2008[ Gardner & Ringhofer 


199G 


Haas 2011 Manfredi & Haas[ |2001| Manfredi 


2005 


Melrose[ 20081 Mendonca 2011 Shaikh & Shukla 


2007 


Shukla & Eliassonl 2006 ( 2007 2010 ,2011 Tsintsadze & 


Tsintsadze[ |2009|, |Vladimirov & Tyshetskiy |20TT 


Wil- 


helm 1971 ) associated with the quantum statistical elec- 



tron pressur e ( Chandrasekhar 1931 1939 Landau fc Lif- 
shitz| |1980 ), the quantum Bo hm potential (jGardner'fe 
Ringhoferl 11996] [Haas[ |20TT| |Manfredi fc Haas[ |2001[ 
Manfredi! [20051 IShukia & Eliassonl 120101 120111 IWii- 



helm[ 1971 ) through which degenerate electrons can tun- 
nel through (often known as the quantum electron recoil 
effect), as well as the electron-exchange and electron- 



correlations 


potentials ( 


Brey et al. 1990 


et al.\ 2008 




Hedin & Lundqvist 


1971 


)■ 



Crouseilles 



It has been 

found that the above mentioned quantum forces acting 
on degenerate electrons in quantum plasmas introduce 
new dispersive features to electron plasma oscillations 



(Bohm 1952! Bohm & Pines 1953 Klimontovich & Silin 



1952a|b ) with frequencies in the x-ray regime, which can 
be assessed by using collective x-ray spectroscopic scat- 



ter i ng tec hniques (| Glenzer fc Redmer 
al 



20071. In fact, Glenzer et al. (2007) have reported 



2009 Glenzer et 



observations of electron plasma oscillations (EPOs) in 
warm dense matter (with the peak electron number den- 
sity r~j 3 x 10 23 cm -3 and the equilibrium electron and 
ion temperatures of 12 eV (~ 1.5 x 10 5 degrees Kelvin). 
The latter is different from the electron Fermi tempera- 
ture T F = {h 2 /2k B m e )(3Tr 2 n ) 2/3 1.7 x 10 5 degrees 
Kelvin, corresponding to an electron number density 
no ~ 2.5 x 10 23 cm -3 ). Thus, the recent experiments of 



Glenzer et al. (2007) have unambiguously demonstrated 
the importance of the quantum statistical pressure and 
quantum electron recoil effects on the frequency spectra 



of EPOs ( 


Bohm 


1952 


Bohm & Pines 


1953 


Klimon- 


tovich & S 


ilin 


1952a b 


Shukla & Eliasson 


2010[ 20111, 


although a previous experimental investigation Watan- 



abe ( 1956 ) had already established the quantum disper- 
sion properties of EPOs in metals. 

Very recently, |Shukla fc Eliasson] ( 2012a|b|c ) discov- 
ered an oscillating shielded Coulomb (OSC) potential, 
(also referred to as the Shukla-Eliasson attractive po- 
tential (SEAP) (jAkbari-Moghanjoughi[ |2012[)), which is 



valid for the potential energy much smaller than ksT 
and mc 2 , around a stationary test ion in an unmagne- 
tized quantum plasma, where c is the speed of light in 
vacuum. The SEAP arises due to collective interactions 
between an ensemble of degenerate electrons that shield 
an isolated ion at atomic dimensions. The profile of the 
SEAP in quantum plasmas resembles the Lennard- Jones 
(LJ) potential in atomic gases. The Shukla-Eliasson at- 
tractive force (SEAF), defined as minus the gradient of 
the SEAP, brings ions closer to form ion clusters in quan- 
tum plasmas. In this paper, we demonstrate the forma- 
tion of ion clusters at atomic dimensions by perform- 
ing computer simulations of the equations of motion for 
an ensemble of ions that are interacting with each other 
through the SEAF. 



I. THE SEAF 

The existence of the SEAP, which is obtained from 
Fourier transformation of Poisson's equation with the 
quasi-stationary electron density perturbation deduced 
from the linearize d continuity and generalized mo- 
mentum equation ( |Shukla fc Eliasson 2012a|b|c ) for 
non-relativistic, degenerate electrons in a dense quan- 
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FIG. 1: The value of a as a function of as/ro. The critical 
value a = 1/4 is indicated with a dotted line. 




a=0.627 

a=0.25 

a=0 



FIG. 2: (Color online) The electric potential as a function 
of r for a — 0.627 (solid curve), a = 0.25 (dashed curve) and 
a — (dash-dotted curve). The value 0.627 is the maximum 
possible value of a in our model, obtained for as/ro « 0.15. 



turn plasma, critically depends on the electron number 
density through the parameter a — h 2 uj 2 e / Am 2 u^, 



where 



(</3 



1/2 



The parameter a mea- 



sures the quantum electron recoil effect caused by 



the quantum Bohm potential (Gardner & Ringhofer 
19961 IManfredi fc Haas[ |2001[ |Manfredi[ |2005[ |Wff 
helm] |l971[ ) V B = (7i z /2ro)(l/- N /u)V^y / re compared to 



the quantum statistical Fermi electron pressure and 
the electron-exchange and electron-correlation effects 
arising from 1/2-spin of the degenerate electrons. Here 
f* = hfin 2 ) 1 / 3 /mr Q is the electron Fermi speed and 

v ex = (0.328e 2 /mr ) 1/2 [l + 0.62/(1 + 18.36a B nl /3 )}^ 2 
includes the effects of electron exchange and elec- 
tron correlations, where tq — n represents the 
average inter-electron distance. The expression for 
v ex is derived by linearizing the sum of the elec- 
tron exchange and electron correlation potentials 



( Brey et al. 
0.985e^n i / J [l 



1990 



Hedin fc Lundqvist[ |1971| ) V xc = 
{0.034/a B n 1 / a )ln(l + IS. 2,7 a bU 1 ' 3 )} . 
that a depends only on as/ro, as 

0.62/(1 



We note that a depends only on 
a ~ 9.37i Wr )/fl + (37r 2 ) 2 / 3 (q B /r ) 
18.36a B /r )] 2 . [Shukla fc Eliassonj p012a|b|c|) found that 



attractive potentials between ions exist only for a > 1/4 
Figure 1 displays the value of a as a function of as/ro, 
where one observes that it is above a critical value 
0.25 only for a limited range 2 x 10~ 2 < a B /r < 1, 
corresponding to an electron number density in the 
range 5.4 x 10 19 cm~ 3 < n < 6.7 x 10 24 cm~ 3 (with 
as = 5.3 x 10~ 9 cm). The maximum value is a ~ 0.627 
at ag/r 0.15, corresponding to the electron number 
density no ~ 2 x 10 22 cm~ 3 , a few times below solid 
densities. The validity of the SEAP has been further 



expanded ( Akbari-Moghanjoughi 2012) for wider den- 



sity ranges by including Chandrasekhar's (1931 1939) 



generalized pres sure l aw for degenerate electron fluids 
and Salpeter's (1961) electron-exchange and electron- 



correlation potentials that are of astrophysical interest 
(e.g. the cores of white dwarf stars). The nonlinear 
shielding effects on the SEAF is discussed by |Shukla ei 
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FIG. 3: (Color online) The initial positions of ions in the 
particle simulations. 

For a > 0.25, the profile of the electric potential as a 
function of distance r around a stationary test ion charge 
Q is ( |Shukla fc Eliasson[ |20T2a|b|c| ) 



4>{r) = — [cos(kir) + sin(fcir)] exp(— k r r), (1) 

which is referred to as the SEAP. Here Q is the ion 
charge, 6* = 1/V4a - 1, h = (k s /y/4a)(y/4a - 1) 1/2 , 
and k r = (k s /y /r Aa)(y /r 4a + l) 1 / 2 , with k s — uipe/u* be- 



ing the modified inverse TF screening length. The spatial 
profile of the SEAP in Fig. 2 shows a distinct minimum 
for the case a = 0.627. The negative part of the SEAP, 
given by Eq. (1), resembles the LJ potential and leads to 
a short-range SEAF between neighboring ions. On the 
other hand, the po tential distribution around a t est ion 
for a < 0.25 reads ( |Shukla fc Eliasson| |2012a|b|c[ ) 



(r) 



Q 

2r 



[(1 + b) exp(-fc+r) + (1 - 6) exp(-k-r)], (2) 
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FIG. 4: (Color online) The positions of ions t — 15000 for a) a = 0.627 and b) a = 0.25, showing the clustering and solidification 
of ions for a — 0.627, but not for a = 0.25. 



where b = — 4a, and k± = fc s (l =p 

y/1 - 4:a)^ 2 /V2a. For this case, the potential is posi- 
tive and monotonically decreasing (cf. Fig. 2), giving 
rise to only a repulsive force (similar to the TF force) 
between ions. In the high- or low-density limit, where 
a — > 0, we recover the modified TF screened Coulomb 
potential </>(r) = (Q/r) exp(— k s r). 

II. DEMONSTRATION OF ION CLUSTERING 

We here present a computer simulation study of the 
dynamics of a system of ions interacting with each other 
under the action of the SEAF. For this purpose, we nu- 
merically solve the equations of motion for a system of 
ions with equal charges and masses, given by 

M ^ = v ^(l R d) - VjVcfo) - Mw jt (3) 

where instantaneous position of each ion is determined 
from chcj/dt = Vj. Here FLy = — Tj is the ra- 
dius vector between particle i and j, rj(t) the position 
and Vj(t) the velocity of the jth ion, M the ion mass, 
Vj denotes the gradient of cf) at position r J; and v de- 
notes an effective collision frequency, which tends to re- 
tard the ion motion. The external confining potential, 
V c (r) — (M /2){w']_r']_ + tu 2 z 2 ), of charged particles may 
have different amplitudes wj_ and uj z perpendicular and 
parallel to the z-axis, respectively, where r 2 ^ = r 2 + r 2 . 

In order to demonstrate the clustering of ions un- 
der the SEAF, we now carry out particle simulations of 
Eq. (3) with 1000 particles, initially randomly placed 
in space, as shown in Fig. 3. In the first set of sim- 
ulations, displayed in Fig. 4, we consider interactions 
of the ions in the absence of the external confining po- 
tential V c , viz. ui z = u)± = 0, and with v = 0.01 t -1 . 



The positions of ions at the end of the simulations are 
shown in Fig. 4 at time t = 15000 r. Here the po- 
sitions and time are in units of L s = k^ 1 and r = 

Af 1 /2/Q fc 3/2 = ( 7r /l 6 )l/8( OB / ro )3/8 w -l /a 3/8 Zi) respec _ 

tively, where w P i — (AnC^niQ/M) 1 / 2 is the ion plasma 
frequency, and n^o the equilibrium ion number density, 
related to the electron number density no via the quasi- 
neutrality condition ZiUiQ = no, where Zi is the ion 
charge state. For a = 0.627, which leads to a poten- 
tial minimum (cf. Fig. 2), we observe the clustering 
of ions and the formation of large-scale ionic structures. 
The clustering of ions is a relatively slow process in com- 
parison with the ion plasma period 2ir/ujp i . Ion pairs and 
smaller clusters are initially formed, and later the larger- 
ion clusters are gradually formed by the agglomeration of 
smaller ion clusters. As a contrast, for a = 0.25, we see in 
Fig. 4 that there does not exist condensation/coalescence 
of ions. This is due to the fact that there is no potential 
minimum for this value of a (cf. Fig. 2), and hence the 
inter-ion force is always repulsive for this case. 

Furthermore, we have carried out a set of simula- 
tions with a symmetric potential V c with lo±_ = u> z = 
2 x 10 _3 t _1 . Here, as seen in Fig. 5, almost spher- 
ical non-Coulombic ion crystals are formed for both 
a = 0.627 and a — 0.25. We also performed simula- 
tions including an asymmetric external potential V c with 
cj± = 6 x 10~ 3 r _1 and uj z = 2 x 10~ 3 r -1 , providing a 
stronger confinement in the perpendicular direction. Fig- 
ure 6 displays the final state for a — 0.627 and a = 0.25, 
where ions, in both cases, form non-Coulombic ion crys- 
tals elongated along the z-direction. In general, the for- 
mation of non-Coulombic ion crystals is attributed to 
the balance between the external and inter-ion poten- 
tials, where the system tends to a configuration of a min- 
imum potential energy. We note that similar configu- 
rations have been previously reported for both charged 
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macro-parti cles (|Wuerker et al. 1959) and ions (Drewsen 
Wed] [19981 |M0lhave fc Drewsen| |2000l |Staanum et al.\ 
2010) confined by external potentials in the Paul trap. 
Finally, in Figs. 7 and 8, we continued the simulations 
in Figs. 5 and 6 and set the external confining poten- 
tial to zero at time t = 15000r. For both cases, with 
symmetric and asymmetric potentials, we see that ion 
crystals remained tightly packed for a = 0.627, where 
it performed damped oscillations as elastic solids before 
settling down to the final states in Figs. 7(a) and 8(a). In 
contrast, for a = 0.25, the ion cloud expanded to larger 
and less dense ion clouds, as seen in Figs 7(b) and 8(b) 
at t = 20000r. Hence the short-range potential around 
ions for a — 0.627 (cf. Fig. 2) lead to the solidification 
of the ions, which is not the case for the the short-range 



repulsive potential for a = 0.25. 



III. SUMMARY AND CONCLUSIONS 

In this paper, we have carried out particle simulations 
to demonstrate clustering of ions due to the newly found 
SEAF arising from collective interactions between an en- 
semble of degenerate electrons that shield ions in dense 
quantum plasmas. Specifically, the SEAF leads to clus- 
tering/condensation or coagulation of ions in the absence 
of an external confining potential for charged particles. 
However, ion clustering can be put on the firm footing 
by calculating the dynamical ion structure factor (DISF) 
based on the fluctuation-dissipation theorem and the di- 
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FIG. 7: (Color online) The positions of ions at t — 20000r for a) a = 0.627 and b) a = 0.25 for the case of a symmetric 
parabolic potential in Fig. 5, where the potential is set to zero at t — 15000r. The ion cluster remains tightly packed for 
a = 0.627, while the ion cloud expands for a = 0.25. 




FIG. 8: (Color online) The positions of ions at t = 20000r for a) a = 0.627 and b) a = 0.25 for the case of an asymmetric 
parabolic potential in Fig. 6, where the potential is set to zero at t — 15000r. The ion cluster remains tightly packed for 
a = 0.627, while the ion cloud expands for a = 0.25. 



electric constant of degenerate electrons and strongly cor- 
related ions in a viscous quantum plasma. It may well 
turn out that the DISF will reveal long-range correlations 
between ions. It is our belief that the formation of ion 
nano-clusters is going to play a valuable role in the area 
of high-density compressed plasmas with degenerate elec- 



trons (|Glenzer & Redmer 


2009 JMalkin et al.\ 2007 


Son 


&Fisch 2004 2005[ 2006|) 


for ICF to succeed, and a 


so m 



the emerging field of nano-material sciences (e.g. nanodi- 
odes, metallic nanostructures for thin films ( |Crouseilles et 



at. 2008 1 , nanowires) , where closely-packed ions will lend 
support to enhanced fusion probabilities (with anoma- 
lous fusion cross-sections) for controlled thermonuclear 



ICF, and may also influence the electric properties (e.g. 
resistivity) of new high density plasma materials at rel- 
atively high temperatures. Finally, we stress that the 
Cooper pairing of ions at atomic dimensions shall pro- 
vide possibility of novel superconducting plasma-based 
nanotechnology, since the electron transport in nanos- 
tructures would be rapid due to shortened distances be- 
tween ions in the presence of the novel SEAF. 
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